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Abstract 


This report summarizes the progress of the research supported by NASA Grant 
NAG-1-643. With the goal of tailorability in mind, the in-plane stiffness characteristics 
of a particular grid stiffened plate configuration under axial and shear loads have 
been studied. The contribution of the skin to the stiffener network and the resultant 
skin/rib interaction is analyzed. For the given plate geometry and loads, it is shown 
that an optimum configuration does exist. 

In order to achieve optimally designed practical plate configurations, buckling 
constraints need to be included in the design. Due to the complex geometry and 
loading of the plates, a simplified local buckling analysis of isolated stiffeners and 
triangular skin elements between the stiffeners is considered. Development of a 
stiffener buckling analysis that represent stiffeners as shear deformable plate 
elements is presented. 


Page 1 


Introduction 


The potential advantages of applying grid stiffening concepts to composite structural 
components are numerous and varied. They include high stiffness-to-weight ratio, 
tailorability, damage tolerance, and economy of manufacture. 

One of the earliest grid stiffening concepts was isogrid, so named for its 
pseudo-isotropic nature. While this may be appropriate for some applications, most 
structures have higher stiffness and strength requirements in specific directions. 
Thus it is expected that the efficiency of stiffened composite plates can be improved 
by tailoring their properties and geometries to meet the requirements of a specific 
application. 

The work being performed for designing minimum weight grid stiffened plates under 
combined in-plane loadings is discussed in this report. The initial effort has 
concentrated on developing appropriate modelling of a basic element of this plate 
configuration in order to understand the response of these plates under the individual 
load cases. A parallel work concentrated on the development of a simplified buckling 
analysis of plate sub-components in order to incorporate practical stability constraints 
into the design procedure. 


Stiffness Characterization 

The work to date has focused on the in-plane stiffness characteristics of a grid 

stiffened composite plate under tension and shear loads. In general, a typical panel 

would consist of a laminated composite skin with a network of attached rectangular 

cross-section ribs, Figure 1. As shown in the figure, a single stiffened cell was 

isolated for analysis. The reasons for this are two-fold. First, in modelling a large 

panel with a repetitive stiffener pattern, geometric symmetry can be exploited to 
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reduce the problem to one of analyzing a single representative cell. Thus, a fine 
mesh can be used to model the component without incurring excessive 
computational expense. Secondly, some thought has been given to examining low 
width-to-length ratio structures, such as beams or ribs, where there may only be a 
single row of cells in one direction. With appropriate boundary conditions the same 
cell mesh can thus be used to model either a large panel or a slender member. 

The particular cell under investigation was 254 mm (10 in) long by 147 mm (5.77 in) 
wide with a 60 degree included angle between the stiffeners in the axial direction. 
The skin thickness and stiffener width remained constant at 2.54 mm (0.1 in) and 19 
mm (0.75 in) respectively. In line with current filament winding manufacturing 
techniques, the ribs were composed of unidirectional material. The skin was 
quasi-isotropic, with a (0/ + 45/ -45/90) layup. Elastic properties were representative 
of a typical graphite-epoxy material. 

Since the objective of this task was to study rib/skin interaction, a design variable 
'stiffener ratio' (SR) was defined to be the ratio of stiffener height to skin thickness. 
The stiffener ratio of SR = 0 represents an unstiffened rectangular plate, while SR -* 
oo approaches a lattice grid (no attached skin). In practice, the lattice grid was 
modelled by setting the skin properties equal to zero. 

Axial Loading Results 

In this portion of the study, it was desired to model a representative portion of a large 

grid stiffened panel under uniform axial tensile load. As such, particular attention 

was paid to achieving the proper boundary conditions. To apply the axial load, a unit 

displacement in the x direction was imposed on all the nodes at one end while 

leaving them free in the y direction. However, it was found that the longitudinal 

edges had a tendency to curve or bow outward, see Figure 2. Caused by the ribs 

attempting to pivot or 'scissor' about the central node, the effect increased in severity 
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with larger stiffener ratios. While this behavior may be representative of a cell with 
free edges, for an internal section from a panel it violates displacement compatibility 
between neighboring cells. To eliminate the curvature, nodes along the longitudinal 
edges could simply be fixed in the transverse direction but left free in the axial 
direction. However, this restraint is probably not indicative of most applications 
(unless the panel were contained within a rigid frame) and results in artificially high 
axial stiffness values. 

To model the longitudinal edges properly, it was desired to impose the condition that 
they remain straight, but still allow contraction/expansion in the transverse direction. 
This was achieved by adding a row of elements to each side of the plate, see Figure 
3, with high transverse (E 2 ) and shear (G 12 ) moduli to reduce the edge curvature, but 
a very low axial (£,) modulus so as to not alter the axial response of the panel. It was 
found that imposing the straight edge constraint with these beam like features 
increased the axial stiffness of all stiffened panels. The effect was more pronounced 
for higher stiffener ratios, rising to a maximum 5.4% improvement for SR = 10 (Figure 
4). 

As shown in Figure 4, the axial stiffness increased with thicker stiffeners, as 
expected. However, it is the structural efficiency on a weight basis, or equivalently 
the stiffness per unit volume (specific stiffness) which is of interest. Except for very 
low stiffener ratios, the efficiency of the panels decreases appreciably with increasing 
stiffener ratio, as indicated in Figure 5. At a stiffener ratio of SR = 10 the specific 
stiffness is only 50% of that for an unstiffened panel. Consistent with this trend, the 
specific stiffness for a lattice grid was extremely low at 0.0364 N/mm 4 , which 
represents just 4.4% of the unstiffened plate value. 

The panels appear to be relatively flexible under axial loading because of the 

stiffeners' tendency to scissor together. Therefore, as the stiffener ratio becomes 

larger, the skin is less effective in controlling this action. Transverse contraction of 
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the panel became increasingly pronounced with the higher stiffener ratios, and in 
particular the lattice panel (Figure 6). 

At very small stiffener ratios the specific stiffness increased with increasing stiffener 
ratio and reached a peak, indicating an optimum value of the stiffener ratio. It is 
believed that the small increase in specific stiffness at the low stiffener ratio was due 
more in part to laminate ply synergism than to geometric stiffening effects. Figure 7 
shows that the addition of a small amount of 30 degree material to the quasi-isotropic 
substrate results in an increased E x for the stiffener material. Therefore, the material 
itself became stiffer, resulting in improved panel stiffness. Unfortunately, this effect 
only appeared at very low stiffener ratios, with the trend to lower moduli giving 
reduced specific stiffness for the panels above SR = 0.25. 

From these results, it would appear that this grid configuration is not well suited to 
axial loading. Since the primary mode of deformation appeared to involve rotation 
of the stiffeners, some means of controlling this would likely show substantial 
efficiency improvements. This might take the form of transverse stiffeners at the 
nodes, which results in isogrid, or employing an anisotropic skin with increased 
transverse stiffness. 

Shear Loading Results 

The problems associated with modelling a panel under shear loads were similar to 
those encountered for the axial case. That is, it was desired to constrain all edges 
to remain straight but still allow free contraction or expansion in both the axial and 
transverse directions. This was eventually accomplished by adding a second pair of 
rows of elements to the ends of the panel. The shear loading was applied by 
imposing a unit displacement in the direction of the diagonal, with the opposing node 
fixed (Figure 8). 
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Qualitatively, the results were somewhat similar to those for axial loading. With 
increasing stiffener ratio, the shear stiffness rose concurrently (Figure 9). However, 
the peak shear stiffness for SR = 10 was 160% higher than that for SR = 0, whereas 
the equivalent increase for the axial case was only 90%. On a weight (volume) basis, 
the specific shear stiffness again decreased significantly after an early peak (Figure 
10). Although the improvement at low stiffener ratio was considerably more 
pronounced than for the axial case, the increase was still limited to only about 4% 
over the unstiffened plate, at a SR = 0.4. For higher stiffener ratios, the degradation 
in specific shear stiffness was not as severe as that for the axial case. For example, 
at a stiffener ratio of SR = 10 the shear panel retained 69% of the unstiffened plate 
stiffness, versus only 50% for the axial case. The lattice panel had a specific shear 
stiffness equal to 0.155 N/mm 4 , which represents 44% of the plate value. This 
displays a considerable improvement over the lattice grid under axial loads, which 
only retained 4.4% of the corresponding unstiffened plate stiffness. Also, the 
deformed lattice grid does not appear significantly different than the stiffened panels 
(Figure 11), which again differs from the behavior under axial loading (Figure 6). 

In summary, the stiffened panel behaved similarly to the axial panel in a qualitative 
sense. Quantitatively, the shear panel appeared to be more sensitive to the changes 
at low stiffener ratios, but less prone to degradation in specific stiffness with 
increasing stiffener ratio. The lattice grid in particular displayed much improved 
performance over its axial counterpart. 

Summary of Stiffness Characterization 

Imposing suitable boundary conditions on the stiffened cells (i.e. the straight edge 
constraint) was found to have a significant effect on the results obtained. In the shear 
case, the addition of beams was essential to distribute the nodal point loads, 
eliminating gross local deformations. These results highlight the need to model 
specimen boundary conditions as accurately as possible. 
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For both the axial and shear loaded plates, an optimum specific stiffness existed at 
very low stiffener ratios. Although this effect was more pronounced for the shear 
case, the increase was still relatively small. Aside from this modest initial increase, 
the general trend was to decreasing specific stiffness for increasing stiffener ratio. 
The limiting case was provided by the lattice grids, which displayed poor specific 
stiffness, especially so for the axially loaded example. This result suggested that in 
order to increase the efficiency of the grid stiffened panel under axial and shear 
loads, it is essential to control the scissoring of the stiffeners. 

For the rib plus attached skin panel, preventing the scissoring action might include 
using an anisotropic skin with increased stiffness in the transverse direction. In 
addition, the results for very low stiffener ratios hinted that it may be possible to 
improve the stiffener material properties via judicious choices of rib and skin ply 
combinations. For either the skinned or the lattice panel, the addition of transverse 
stiffeners at the nodes should also result in substantial improvement in panel 
stiffness due to triangulation of the structure. Of course, adding transverse members 
of the same thickness as the original ribs would result in isogrid, with its attendant 
pseudo-isotropic behavior. However, if the panel were to support loads of a 
directional nature, it would likely be possible to adjust the thickness or width of the 
ribs individually to arrive at a reduced weight configuration. 


Page 7 



Stability Constraints 


The in-plane stiffness analyses of the previous section were conducted without any 
type of buckling constraints imposed on the designs. Thus, for maximum specific 
stiffness, a design was similar to that of a flat plate with rather low stiffeners. 
However, intuitively one expects that in practice the performance of this structure 
would be poor due to its low out-of-plane bending stiffness, and thus attendant low 
resistance to buckling under compressive loads. Besides the constraints on specific 
stiffness and local failures associated with the material strength failure in the 
stiffeners and the skin panels, constraints against stability failures must be included 
in the design in order to achieve meaningful structures. Global buckling of the panel, 
as well as the local buckling of the skin between the stiffeners and the stiffeners 
themselves must be included in the optimization process. 

In evaluating the suitability of various proposed buckling schemes, it should be noted 
that verification of the model appears to be difficult. Experimental results in the 
literature for buckling tests on grid stiffened composite panels are very limited. The 
data that is reported often exhibits a considerable amount of scatter, generally 
thought to be a result of manufacturing imperfections. Due to the complex geometry 
of the panels, an exact analytical solution does not exist. The best estimate of panel 
response may be that provided by a detailed finite element (FE) analysis, which may 
result in a prohibitively expensive design optimization run. Rather than performing 
a sophisticated finite elements buckling analysis of the grid stiffened plate, the 
approach used in the present study is to use simplified buckling analysis of the 
overall plate and the individual components. Following is the development of the 
procedure for the buckling analysis of the stiffeners. 
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Stiffener Buckling Failure 


Grid stiffened fiber-reinforced composite plates are primarily manufactured as 
filament-wound structures. The filament winding process as an economical method 
of fabrication makes the grid stiffened structural concept attractive. For stiffeners 
with unidirectional fibers, the stiffener material is orthotropic with respect to the 
stiffener axis. For typical high performance composite materials, the shear stiffness 
of such a configuration will be low compared with the longitudinal modulus of the 
material. Even for stiffeners which are laminated with plies of various angular 
orientations, the transverse shear stiffness is low. The low transverse shear 
stiffness, combined with the possibility of low ratios of stiffener length to thickness 
ratio, point to the importance of the effects of transverse shear deformations in 
computing the buckling load for a stiffener. With this point in mind, a procedure is 
developed to predict the critical end load for a composite stiffener. The stiffener is 
assumed to be a simple blade type, and thus is treated as a rectangular plate with 
one edge joining the skin, the opposite edge away from the skin free, and the other 
two edges loaded uniformly along the stiffener axis. The loads for the stiffeners are 
calculated from the finite analysis of the plate already performed for the stresses and 
displacements. The boundary conditions at the loaded edges of the stiffeners are 
probably between that of simply supported and fixed conditions, with indications in 
the literature that it is closer to a fixed condition. In considering anticipated buckling 
modes, if the stiffeners are assumed to buckle simultaneously with the triangular skin 
sub-elements, then the hub of the stiffeners would undergo no rigid-body rotation, so 
that the stiffener end would be fixed except for a small amount of flexibility in the hub. 
One possible treatment is to assume that the stiffener ends are fixed, but that the 
effective stiffener length is greater than the exposed stiffener length by some fraction 
of the hub diameter. Because many of the earlier works have considered simply 
supported ends than fixed ends, the computational procedure developed here for the 

simply supported end conditions. This will result in an conservative design results. 
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In the following analysis the stiffener plate is considered to be a laminated composite 
material which is orthotropic with respect to the stiffener axis. First order shear 
deformation theory is used in a Rayleigh-Ritz procedure to solve for the critical 
buckling load and the buckling mode shape. 

Buckling Calculations 


Theory 

The stiffener is treated as a rectangular plate supported on three edges, with a 
uniform compressive load, N', applied at the two opposing supported edges. The 
plate dimensions and coordinate axes are shown in Figure 12. The plate material is 
assumed to an orthotropic composite laminate. The plate theory used is the 
first-order transverse shear deformation theory for plates, in which shear deformation 
is assumed to be linear through the thicknesses. 

The assumed displacement field in the plate is given by 

u,(x,y,z) = u(x,y) + Z</> x (x,y) 

u 2 (x,y,z) = v(x,y) + Z<p y {x,y) (1) 

u 3 (x,y,z) = w(x.y) 

where u u u 2 , and u 3 are the displacement components in the x, y, and z directions, 
respectively, u, v, and w are the corresponding displacements of the mid-plane of the 
plate, and <f> k and <f> y are the rotations due to bending of a material line initially 
perpendicular to the mid-plane, about the y- and x-axes, respectively. (Using the right 

A 

hand rule, the sense of the rotations are j for cf> x , and — / for <f> y .) The 
strain-displacement relation for small displacements, 

£ ij = \ i U U + U 1 ,/) 


( 2 ) 
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results in the following expressions for strains: 



(3) 


Energy considerations can be used to derive the differential equations of equilibrium 
for the plate. Assuming a linearly elastic material, the strain energy of the plate, U, 
can be written as 


U = 


. „ h 

Jf 

J Q J -h 


—o^dzdQ 


( 4 ) 


where £2 is the domain of the plate in the x-y plane. The potential energy of the 
applied loads, V, is 


V = - f t)UjdS 2 (5) 

*'s 2 

where S 2 is the portion of the plate surface where a traction t is applied. The total 
potential energy, n, is 


n = U + V (6) 

and by the principle of minimum total potential energy, the variation of the total 
potential energy is zero, or 
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an = o 


(7) 


When the loading applied to the plate is limited to in-plane edge loads, the resulting 
system of equations corresponding to virtual displacements of the five displacement 
variables, u, v, w, </>„ and 4> y , are: 


dNy dNy 
dx dy 


dN 6 | dN 2 
dx + dy 


= 0 


dQ i 

dx 


dQ o 

+ -1- + N = 0 

dy 


- dM r 1 

dx 

ax 


a/w„ 


dy 


+ 0 1 = 0 


dM 2 

dy 


+ Q 2 = 0 


( 8 ) 


where the in-plane edge force resultants, /V„ the moment resultants, M it and the edge 
shear force resultants, Q , , are given by 


_ 


rf 


N, 



'i * 




_z _ 



-h 



2 


a,dz 



i= 1,2,6 


and the transverse force term, N, is 


N = N y 


£w 

dx 2 


+ n 2 


a 2 w 

dy 2 


+ za/ 6 


d 2 w 

dxdy 


(9) 


( 10 ) 
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which arises by considering the projection of the in-plane force resultants on the 
z-axis for the case where the plate is deflected. 


To treat the buckling of a plate stiffener, consider the following boundary load 
conditions: 


x = 0, i\ A/ 1 = —N' x N 6 = 0 

y = 0,b: N 2 = 0 N 6 = 0 

From the first two of equations (8), we have the solution that 


( 11 ) 


Ni(x,y) = — N' x = Const 

A/ 2 (x,y) = 0 in Q (12) 

W 6 (x,y) = 0 


so that we can now limit our consideration to the following three equilibrium 
equations: 


£Qi dQ 2 
dx dy 



- dM 1 

dx 


dMe 

dy 


+ (?i = 0 


-dMe 

dx 


dM 2 

dy 


+ Q 2 = 0 


(13) 


Consider the constitute equations for an orthotropic lamina with respect to its 
principle axes. For the transverse shear, we have the relations 


— <*23 — 20^4 — 2G 23 £ 4 

( 14 ) 

<*5 = <*13 — 2C 55 £ 5 = 2G 13 £ 5 

where the symbols will be defined presently. The constitutive equations relating the 

in-plane stresses and strains are those for an orthotropic plate where the transverse 

normal stress is zero, and in summary we have the following constitutive relations: 
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Q-i i O12 0 0 

Ql2 O 22 0 0 


0 



0 

0 

0 


0 

0 2C44 0 0 

0 0 2C55 0 

0 0 0 2Q ( 



( 15 ) 


where the same contracted notation is used for stresses, <7, as was introduced for 
strains in equations (3), and where 


Q11 — 


1 - V12V2/ 


Q22 ~ 


Vl2^2 


i — yi2V2i 


C44 — g 23 


C 55 — ^13 


W12 — . 

1 ~ V12V21 
^66 = ^12 


(16) 


where £, and £, are Young's moduli, y 12 and y 21 are Poisson's ratios, and G a , G 13 and 
G„ are shear moduli, the subscripts having the usual meanings. 


To transform the constitutive relations for the lamina to an arbitrary coordinate 
transformation, 



or 


{X,} = [g {X'y} (18) 

where the primed axes are the principle orthotropic axes. The the transformation for 
a fourth order tensor is 


Tjjkl ^im^jn^kr^ls^mnrs ( 19 ) 

where T ijkl is in the unprimed axis system. Noting that in expanded notation 
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Ol6 — ^1112 > 


etc. 


(20) 


C44-C 


2323 > 


the transformation results in 


<3^ = 0^ cos 4 # 4- z(C? 12 4- 2Q 66 ) sin 2 # cos 2 # 4- C? 22 sin 4 # 

Oi2 ” (^ii 4" O22 — 4Qgg) sin # cos # 4" Qi 2 ( sin 0 4~ cos 0) 

<?2 2 = Qn sin 4 # 4- z(0 12 4- 2Q 66 ) sin 2 # cos 2 # 4- Q 2 2 cos 4 # 

Oi6 = (Oil - 0i2 - 2 c ? 66) sin 0 cos 3 # 4- (0 12 - Q 22 4- 20 66 ) sin 3 # cos 6 

Q 2 g = (On Oi 2 2066) sin 9 cos 9 4” (0i 2 — 0 22 4- 2Qg 6 ) sin 9 cos 9 (21) 

©66 = (On 4" 0 22 — 20 12 — 2Qgg) sin 9 cos 9 4" 066( sin 6 4” cos #) 

Qu = C 44 cos 2 9 4- C 55 sin 2 # 

Ogg = C 55 cos 9 4- C 44 sin 9 
C 45 = (C 44 - C 55 ) cos 0 sin 9 


Substituting in the expression for strains in terms of the assumed displacements, the 
constitutive relations in the unprimed axis system are 



On Oi 2 2Qi6 0 0 

0 2 i 0 22 2Q 26 0 0 

0i6 0 2 6 2066 0 0 

0 0 0 2C 44 2045 

0 0 0 20,5 2C 55 



u, x 4- z<p x , x 

V,y 4” Z<f>y,y 


4" v, x 4" z(4> x ,y 4" $y>*)J 

1 ( i , dw} 

~ [ * y ~: 

± , dw 

^ , + 17 


( 22 ) 


Consider a laminate made of k orthotropic layers, each oriented at some angle with 
respect to the reference axes. Integrate equation (22) over the thickness to get the 
following constitutive relations between the resultants of equations (9) and the 
assumed displacements: 
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An A 12 A 16 S-n S 12 B 16 



A 2 2 

a 26 

8 12 

8 22 

®26 

• 

A66 

#16 

®: 26 

®66 

• 

• 

D11 

D12 

D 16 

SYM 

• 

• 

CM 

CM 

Q 

0 26 

9 

. 

. 

, 

D 6 6 



= K‘ 


A44 A 45 
A45 a 45 


<t>y + W,y 
<t>X + W -X 


u, x 



( 23 ) 


where K* is the commonly used shear correction factor, assumed here to be ( 5 / 6 ), 
and 


Ajj 

fj 

’ 1 ’ 

N 

1 

1 

:a- 

r 

1 


-if* 

z 

*-I* 

- ^*- 1 ) 

3. 

*=1 V -1 

z 2 . 

/ f — 1 



™ rh 

A (j = £ J cfcfe = £c*(/,* - h k _,) /, 7 = 4, 5 

k = 1 Af=1 


where /»*_, and h„ are the z-coordinates of the lower and upper surfaces, respectively 
(positive z being up), of the k-th lamina with respect to the mid-plane of the plate. 


We now restrict consideration to laminates which are orthotropic with respect to the 
reference axes. These include unidirectional and symmetric cross-ply laminates. For 
these laminates, 8^ = 0 , /j = 1,2,6, and A 16 = A a = D 1S = D w = 0 . For an orthotropic 
material, equations ( 23 ) reduce to 
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1H- + A -Ql 
dx + * 12 dy 

du 

A dv 

*r 4 

A22 17 

d<f) x 

d<t> y 

dx 

+0,2 ir 

d<t> x 

3 * r 

dx 

+ ° 22 a y 





M 6 




( 25 ) 


To employ the Rayleigh-Ritz procedure, begin by recreating the variational 
statements used in deriving the governing differential equations. Multiply the three 
equations (13) by the virtual displacements <5w, 64> x , and 34> y , respectively, and 
integrate over the x-y domain of the plate: 


f* 

1 

so, 

dQ 2 


[dx 

dy 

( 

( -dM 1 

dMe 

I 

"£i 

( dx 

dy 

1 

( ~dM 6 

8M 2 

1 

\ dx 

dy 


* dx 2 


]SwdQ — 0 


+ Q 1 <5 <t> x dQ = 0 


+ (?2 ) S(j)yd£l = 0 


(26) 


Integrate by parts and apply Gauss' divergence theorem to get 
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+ ^ - i£r "‘ irW - I**"* 




( dWx 

V ax 


dd<f> x dS(f> x 
M - I +-^— 

ay 


M 6 + b (j> x Q^jdxdy = J M x S<j> x ds 


db<f), 

dx 


dS 4 > y \ f 

-M e H — + HyQ 2 jdxdy = J M y S<t> y ds 


(27) 


where 


S x =(o,-N' x ^jn x + Q 2 n y 

x ^ I 
M y = M 6 n x + M 2 n y 


(28) 


where n, and n y are the in-plane components of the unit normal vector at the boundary 
of the plate. 

For plates with a combination of simply supported (S.S), fixed, and free edges, the 
RHS terms of equations (27) are identically zero, since on all boundaries either the 
respective displacement components are specified and have a zero variation, or the 
resultant generalized force components are zero. We introduce approximate 
solutions for w, <j) x , and 4> y in the following form: 


a/* % % 

W = 2>,(x,y) (f) x = ^ ' fl/fl x ;(x,y) 0yy ~ ^ 

(=1 /=1 i =1 


(29) 


where ^(x.y), 0 x ,(x,y), and 0 y/ (x,y ) are shape functions which satisfy the geometric 
boundary conditions and which are part of a complete set, and A„ B„ and C, are 
constant coefficients to be determined. The variations which appear in equations (27) 


now become 
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( 30 ) 


5w = 


Arffay) 

/= 1 


S4>x = ^<5B,0 x/ (x,y) 
/=1 


7 

<5</>y = Y^ 5C i e yi( X >y) 
/= 1 


When equations (29) and (30) are substituted into equations (27) and the variation of 
each coefficient seen in equations (30) is taken to be independent of the others, the 
following system of equations results: 


N w 

% 


N A 

*y 

N w 



a + !>” 


+ l>r c j 

= A /' x £M t] Af 

/ = 1,2, 

y=i 

y=i 


7=1 

y=i 



N w 


% 

n, 

*Y 




+ 

S) + 

E< c /-° 

/=1,2,. 


y= i 






N w 


N, 

N 6 

4>y 




+ 

i< 8 / + 

E< c /=° 

/= 1,2,.. 


/=i 


7=1 

7=1 



where 


(31) 
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[* C b dll/, d^j 

: J 

"0 "0 


K U 

K ij 


n » J d\l/i d\l/j d\ Ay\ 

_ o \ 55 dx dx 44 dy dy J 

<=K 2 i =rfK 2 ( A 

,13 _„31_f'f e 

Jq **0 




Ls 13 iy31 

“ii ~ 




*22 - 
Kij - 


j d y dx 

'44“^r y : 

de xi de 


‘a = | | A*|^)dycfx 


3« , <w x , de xj dd xj 7 \ 

D "~eT - IT ' + 1 D - ' IT IT + ) dy 


(32) 


4) 4) 


jv -23 _ „32 _ 

*ij 'v/ “ 


/V /•/> 


D 


de xi dd yj . „ 50 X, 30 yy Y 

1 — + ^66 — : ^ — d 


4 4 

t*e /*b 


dx dy 


Kf = 


b / 

( ddyj ddyj ddyj ddyj „ 

( ° 22 ^ ITT + °66 -t: TT + K A “ e yi 0 yj 


Equations (31) can be reduced to a symbolic form as follows which has the form of 
an eigenvalue problem. 


[[^' 1 ]-W' x [/W]]{/V} = {0} (33) 

For the critical buckling condition the determinant of the matrix in equation (33) must 
be zero, and the resulting eigenvalues, N' x , are the approximate critical buckling 
loads corresponding to various approximate buckling mode shapes. The lowest 
critical value of N\ identified, is the approximate buckling load of interest. 

The required two-dimensional shape functions are chosen here to be built as 

products of pairs of one-dimensional shape functions, one being a function of x, the 

other, a function of y. The following notation for indexing and dimensioning is used: 
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* = 4>i(x,y) = iffOOrffy) , m = 1 ,2,...,N WX , 

i = (m — 1)N wy + n , 


n = 1,2,. ..A/, 


wy » 

N — N N 

"w ,y wx ,Y wy 


0X = 0{(x,y) = l%Xx)g { *\y) , m = 1 ,2 A^ xx , 

i — (m — 1)A/^ xy + n , 


n 1 ,2,...A/^ X y , 
= ty^xx ^xy 


( 34 ) 


4> y = 4>yi(x,y) = fS\x)g^(y) , m = 1 ,2 A/^ yx 

/ = (m — 1)A/^ yy + n , 


n — ^,2,...N ^ p yy , 

N <t>y = N <t>yx N <t>yy 


where the superscripts on f(x) and g(y) indicates the associated variable: w, 0„, or 

4>y ■ 


Stiffener Representation 

In the following paragraphs the boundary conditions for a simply supported loaded 
end is considered, the shape functions are identified, and the resulting expressions 
for the matrix elements (equations (32)) are listed. In those expressions, 3 , is the 
Kronecker delta function, for which 

(1 i —j 

^ 1 (35) 

10 i*j 

and for a typical matrix element, we use the notation 

(36) 

implying that / is a function of m and n by the relations given in equations (39), and j 
is a function of r and s in the same manner. 


The boundary conditions used for this case are 
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w(0,y) = w(«f, y) = w(x,0) = 0 

4> x {x,0) = 0 

0 y (o,y) = 4> v {t, y) = o 


( 37 ) 


and assuming a nominal half-wave buckling mode shape, the additional conditions 
are specified to discount any contributions not symmetric in x about the line x = L/2: 

+ A, yj = - A, y) 

<t>x(j + A, y) - - - A, y) (38) 

+ A > yj = ^y(y ~ yj 

The following shape functions are used: 

•Ai(^.y): f m (x) = sin(2m - 1)^y ; g n (x) = y n 
0*/(*>y): W = cos(2m - 1)*y ; g„(x) = y" ( 39 ) 

S yi (x,y): Ux) = sin(2m-1)jry ; g„(x) = y (n-1) 

and when the integrals of equations (32) are evaluated, the following expressions 
result: 


Ki'=Kmnn = K S A 


55'1 


/i + /< V' 


44'2 


(40 - a) 


where 


/i=(2m-1) 


2 £ b {n+s+:) 

2/ (/7 + S+1) 


£ nsb <B+> - 1) 
2 (n + s-1) 


<5 


mr 




mr > 
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Ki 2 = Kl 2 nrs = K 2 A 55 (2m-1)l--± 


(n+s+1) 


2 (n + s + 1) 


*f 13 i /2 a ^ 

n ‘mnrs ~ A "44 


nb 


(n+s-1) 


2 (n + s - 1) 

^mnrs ~ ^1 1 + ^66^2 = K <^55/3 


where 


/ 1 = [(2m - 1)/r] 


2 1 


(n+s+1) 


2if (n + s + 1) 


f £ nsb (n+s ~ 1) c 

2 2 (n + s- 1) ° mr 

, £ 6 (n+s+1) . 

3 2 (n + s + 1) mr 


K 


23 

mnrs 


— D 12 /i + D 66 / 2 


where 


/, = - (2m - 1> 

< 2 -( 2 n>- 1 )f 


2 (n + s-1) 
(n + s-1) mr 


6 


mr 


K 


33 

mnrs 


— D 22 1*1 + Dgg/ 2 + K 2 A 44 / 3 


where 



£ (n - 1)(s - 1)/7 (n+s ~ 3) 
(/7 + s - 3) 

0 , n = 1 or s = 1 


'mr > 


n > 1 , s > 1 , 


/ 2 = [(2m - 1 )tt] 


2 1 


, (n+s-1) 


2<f (n + s - 1) 


'mr 


i 


.(n+s-1) 


/, = - Z d mr 

3 2 (n + s-1) mr 


(40 - b) 


(40 - c) 


(40 - d) 
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Summary of Buckling Considerations 


The equations governing the buckling of a uniaxially loaded orthotropic stiffener were 
developed using a first order shear deformation theory for plates. An approximate 
solution procedure was developed and implemented on a computer using the 
Rayleigh-Ritz method, for plates with simply supported loaded ends. The method 
also adapts easily to the treatment of uniaxially loaded orthotropic plates with other 
boundary conditions. This method can be used in conjunction with a stability analysis 
of triangular sub-elements of the panel skin to form a local buckling criteria for 
designing grid stiffened plates. Currently, buckling analysis of panel skin 
sub-components between the stiffeners is under way. 
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